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Abstract 

Partition functions of two different matrix models for QCD with chemical potential are computed 
for an arbitrary number of quark and complex conjugate anti-quark flavors. In the large-iV limit 
of weak nonhermiticity complete agreement is found between the two models. This supports the 
universality of such fermionic partition functions, that is of products of characteristic polynomials 
in the complex plane. In the strong nonhermiticity limit agreement is found for an equal number 
of quark and conjugate flavours. For a general flavor content the equality of partition functions 
holds only for small chemical potential. The chiral phase transition is analyzed for an arbitrary 
number of quarks, where the free energy presents a discontinuity of first order at a critical chemical 
potential. In the case of nondegenerate flavors there is first order phase transition for each separate 
mass scale. 
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1 Introduction 



The idea to use random matrix models as a simple model to study the non-perturbative phenomenon 
of chiral symmetry breaking in QCD |T| has been very successful. Apart from direct studies through 
lattice QCD or other effective models, it has become one of the available tools in this area, and we refer 
to [2] for a review on the matrix model approach. In the broken phase at zero temperature and chemical 
potential the applicability of a matrix model has been completely understood by rederiving part of its 
results from the underlying effective chiral Lagrangian picture, that describes the pseudo-Goldstone 
fields. This has been achieved for the partition function the spectral correlation functions of Dirac 
operator eigenvalues [UUHE], as wen as f° r individual eigenvalue correlations [7]. 

A random matrix model that includes the effect of a baryonic chemical potential was introduced 
in [H], studying the nature of the quenched approximation and the chiral phase transition. The same 
model was enlarged to include the effect of temperature, and the phase diagram of QCD with two 
light flavors was predicted [5], including the existence of a tricritical point. The analysis was repeated 
very recently in JOJ) distinguishing between baryon and isospin chemical potential. Only the matrix 
model with isospin chemical potential can be related to an effective chiral Lagrangian of quenched 
QCD [TT] so far. 

The unquenched matrix model partition function with chemical potential for a single flavor has 
been analyzed in great detail, including its analytic solution and the behavior of its zeros [121 113) . 
It has been used as a test case in jlril I14| for lattice algorithms for chemical potential. A complete 
and detailed solution of the unquenched matrix model partition function with several, non degenerate 
flavors has been lacking so far, apart from a first attempt in [E]- It is one of the purposes of this article 
to provide such a solution, both for finite- N as well as in the large- N limit. The characteristic feature 
of the model is the nonhermiticity induced via the chemical potential, which renders the eigenvalues 
of the Dirac operator complex. Under these conditions one has to distinguish in the large- N limit 
between the regimes of weak and strong nonhermiticity (see |17| for a review). 

In ^Hj a random matrix model with complex eigenvalues was proposed for the phase with broken 
chiral symmetry . The resulting spectral correlation functions were computed, and the analytic pre- 
dictions were confirmed by comparing them with the results from quenched lattice simulations |19j . 
both in the limit of weak and strong nonhermiticity. However, the equivalence in the phase with 
broken symmetry between the model proposed in and the original model with chemical potential 
proposed in [H] was only conjectured [THllTo] . 

In a recent paper 0| the authors managed to derive the spectral density of complex eigenvalues 
in the regime of weak nonhermiticity directly from the effective chiral Lagrangian for quenched QCD 
combined with the model [Hj. Exploiting a variant of recently suggested exact replica method |20| 
they arrived to a density profile slightly different from earlier results of ^H]. Both results agree 
asymptotically, and the difference was too small to be distinguished from the lattice data for the 
values of chemical potential used. Under these circumstances it is conceptually important to be able 
to prove the universality of results within the random matrix model approach, apart from matching 
them with first principle lattice data. Without such a universality random matrix models loose much 
of their predictive power, being deduced from global symmetry arguments alone. When the chemical 
potential is absent, universality was proved in [211 122j . For correlation functions involving complex 
eigenvalues only partial results exist for non-chiral random matrix models [23] at weak nonhermiticity. 

One of goals of the present paper is to clarify the issue of equivalence and thus possible universality 
of the two different random matrix models [S] and ^S] for QCD with chemical potential, at the level 
of the corresponding partition functions. We are going to demonstrate that both models agree at the 
regime of weak nonhermiticity for any number of quarks and conjugate anti-quarks, and in that sense 
they are universal. At the regime of strong nonhermiticity the agreement, however, persists only for 
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an equal number of such flavors. For a general flavor content, including only quarks, the two partition 
functions jSJ agree only to the leading order term in an expansion in small chemical potential. We 
will also relate our findings to the so far open question of universal spectral eigenvalue correlations, 
by mapping it to the problem of universality of the so-called bosonic partition functions. 

The object of our investigation can be phrased also in more mathematical terms. In the presence 
of quark flavors the problem of calculation of the random matrix model partition function amounts 
to computing the expectation value of a product of characteristic polynomials (also known as spectral 
determinants) in the model with zero flavors. Characteristic polynomials in Hermitian random matrix 
models have received a lot of attention recently, partly due to their relevance to the behaviour of 
Riemann zeta-function suggested in |24j . On the other hand, such polynomials find important appli- 
cations in theory of disordered and chaotic systems, see [123] and discussion and further references in 
I25J. Different formulas for arbitrary products [22 EFJ E3 and ratios E01 EH E3 EH1 EH ESj of 
characteristic polynomials have been derived and the universality of these expressions has been shown 

HSU2HEZ]. 

Again much less is known for complex matrix models |381 1391 14U| . Here one has to distinguish 
between characteristic polynomials and their complex conjugates. A closed determinant formula for 
arbitrary products of characteristic polynomials and their conjugates (of not necessary the same num- 
ber) for quite a general class of models has been given in |39| . see also related objects emerging in the 
theory of quantum chaotic scattering |41| I17j . Our aim here is thus to compute and compare such 
products within the two models [HI El- These results may also be useful when several sets of replicas 
are needed in the computation of two- or higher fc-point eigenvalue correlation functions, generalizing 

We derive a compact, new expression for the partition function of the matrix model jH] at finite- 
A with arbitrary many quark flavours . It allows us to analyze the chiral phase transition in more 
detail, in particular concerning the influence on nondegenerate quark masses 1 . We find that the first 
order phase transition found in [SI [HI persists, and that it is always driven by the flavor with the 
smallest (or zero) mass. For different mass scales present there is a discontinuity of first order for each 
different mass. This is due to the fact that roughly speaking the partition function can be written as 
a determinant over single flavor partition functions. A similar phenomenon of having two first order 
lines for two flavors was observed in JJJ1; where two different chemical potentials for each flavor were 
introduced. 

The outline of the article is as follows. In the section [2] we define the two matrix model partition 
functions for QCD with chemical potential and compute them for finite- A. We distinguish between 
the presence of only quarks in section l2~Tl and quarks with additional complex conjugate anti-quarks 
in section l2~2l Several technical details of the derivation are summarized in the appendices lAl and iBl 
In section [31 we turn to the large- A limit, where we distinguish between the limit of weak and strong 
nonhermiticity in sections 13. II and l3. 21 The resulting consequences for the universality of matrix model 
partition functions are discussed in section 13.31 In section [I] we exploit the results for finite- A from 
section [21 to investigate the chiral phase transition at a critical chemical potential, for an arbitrary 
number of quarks. In section [SJ we summarize our findings. 

1 We recall that the model is always in the broken phase by definition, it cannot reach the phase transition. 
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2 Matrix model partition functions for finite-TV 



2.1 Partition functions with Nf quark flavors 

We start with defining the first matrix model for QCD with chemical potential, initially introduced 
by Stephanov [Sj . In the sector of topological charge 2 v it is given by 

N f I m/ljv i$ + iA N \ 
-Z^'" 5 (/*;{"»/}) = /^d^JJdet exp -A^g^Tr^ . (2.1) 

/=i Vi^ + Mllv rn f l N+v ) 

Here $ is a complex matrix of size iV x (iV + v). Apart from the unity matrix l n of size n x n we 
have also introduced the rectangular unity matrix of size N x [N + ^) 

{1n) ^ = { j = N + l,...,N + p ■ (2 - 2) 

Eq. (|2.1() contains Nf quark flavors with real masses m/. The chemical potential is added to the 
Dirac matrix of the usual matrix model in a standard way by shifting Ij) — > Ift + jofi. The model has 
the same global symmetries as QCD with gauge group SU (N c > 3) in the fundamental representation. 
The Gaussian weight replacing the gauge field average was chosen for simplicity in [Sj. It will allow 
us to exactly solve the partition functions for finite- N. The variance of the random matrix entries is 
such that the Banks-Casher relation for [i = is satisfied. We will not address the issue of universality 
by allowing for a more general weight function, — » V with V being a polynomial, as in 
|211 122 j . Instead, we compare to a different model given in terms of complex eigenvalues |18| defined 
below, which reduces to the same model P at \i = 0. 

In the presence of /i 7^ a diagonalization necessary to obtain complex eigenvalues of the Dirac 
operator does not any longer amounts to a simple procedure. More precisely, the angular variables 
used in the singular value decomposition of the matrix $ in the form $ = [/lAc^with U\i being 
unitary, no longer decouple in the matrix integral. Although a Schur decomposition to an upper 
triangular form <f> = U(Z + R)U^ remains possible, it does not reveal a natural relevant degrees of 
freedom in the matrix integral. In particular, proceeding in this way one retains the eigenvalues in Z 
complex even after setting \i = 0, which is not a natural choice of integration variables. 

The standard fermionic approach to compute the partition function remains a viable alternative. It 
consists in replacing the determinants by equivalent Grassmann integrals, and further integrating out 
the matrix $ explicitly, and finally performing a Hubbard-Stratonovich transformation. The details 
are given in the appendix and include a more general case of additional conjugate anti-quarks, to 
be addressed below. For real quark masses the result simplifies considerably after choosing a polar 
decomposition, with the Jacobian computed in appendix El After these manipulations we arrive at 
the following result 

Nf 

Jo f=1 

dV J dU det[U^] v exp N(qq) 2 Tr (MUVrV^ + VrV^tf Af) . (2.3) 

Here, dU and dV denote the Haar measure over the unitary group of size Nf x Nf. The matrix 
M = diag(mi, . . . ,niN f ) contains the quark masses and r = diag(n, . . . , rjyy) is the diagonal matrix 

2 Without loss of generality we restrict ourselves to v > throughout the following. 
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Nf 



denotes the Vandermonde determinant. Due to the 



of radial coordinates. A^Jr 2 ) = f] 
unitary invariance we can shift U -> UV* and then rename U — > U* . The resulting unitary integrals 
can be performed exactly using [32], together with the fact that M = M*. Our first main result valid 
for fmite-iV is thus reading 



Z 



(N f ,u) 



N f „ fOO ' 

OxKm,}) ~ e-vWXf^f \[dr k rl + \rl-^) N e- N ^l A Nf (r 2 ) 

J £ -| 



1 



. . det [l„(2iV<^) s 

«j=l,...,A r f 



(2.4) 



Such a compact expression reducing eq. (|2.1j) to Nf real integrations that factorize has previously 
been known only for Nf = 1 ^2]. The integral representation for Nf = 2 in j£] using a Schur 
decomposition of the matrix Q = UVrV^ is more involved. Eq. (|2.4|) is now amenable to a saddle 
point computation, both at weak and strong nonhermiticity. Furthermore, it will be useful when 
investigating the phase transition in section |IJ 



We now turn to an alternative random matrix model for QCD with chemical potential, introduced 
in ^H] in terms of N complex eigenvalues. This type of model can also be solved exactly at any 
finite N. Here we do not have to distinguish between real and complex masses. Moreover, due to 
the powerful technique of orthogonal polynomials the corresponding correlation functions of complex 
eigenvalues also can be found [THj . The model is defined as 



;W» = /II \d 2 z j w(z j ,z])l[m f '(z] + m}) \ 



\&n(z 



2M 2 



w(z, z ) = \z\ T exp 



N 
1^72 



-(z 2 + z 



2*- 



(2.5) 
(2.6) 



The parameter r G [0,1] appearing in the Gaussian weight function w(z,z*) controls the effective 
degree of nonhermiticity. It allows to interpolate between models with real and maximally complex 
eigenvalues for r = 1 and r = 0, respectively. The two partition functions eq. Q2.1j) and eq. (|2.5[) are of 
course different in general at finite N. Only in the Hermitian limit r — ► 1, \i — ► they are exactly the 
same under the identification of the masses: m,fu = rrif j{qq)\/2. In fact, in this paper we will show 
a much more interesting relation: two models are also generally equivalent at an appropriate large N 
limit, under an appropriate correspondence among the relevant parameters (the mapping between the 
two models actually depends on the way the large ./V limit is performed). 

The massive partition function eq. Q2.5|) can be evaluated due to the following observation. We 
can write 

I N N f v 
Z^ u \r; {m f }) ~ ( ft J[ m) (z] + m}) ) , (2.7) 
\j=i f=i I 



where the expectation value is taken with respect to the zero flavor partition function, Zj^ f °' u \t). 
The relation to a product of characteristic polynomials is evident. 

The relevant set of orthogonal polynomials with respect to the weight eq. 
complex generalization of standard Laguerre polynomials, see ^5] : 



2.6(1 are given by the 



p'V) 



1W! [ I) L 



Nz^ 
~27 



(2.f 
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which are given here in the monic normalization. Using the theorem proven in j39j the expectation 
value eq. 1)2.7)1 can be conveniently expressed as a determinant of size Nf x Nf , with entries being 



these orthogonal polynomials. Taking into account the extra factors Yl 



N, 



/=l "7 



m f we arrive 



at [SHI 



r;{m/}) 



det fc ;= 



l,...,N f 



fc-i(-W) 



(2.9) 



We would like to note that the partition function is real for real quark masses as it should be. This is 
not obvious at all when looking at the definition eq. (|2.5j) . The obtained result is exact for finite- N. 
It can also be continued to complex quark masses without modification. 



2.2 Partition functions with quarks and complex conjugate anti- quarks 

In this subsection we enlarge the flavor space by adding pairs of complex conjugate quarks to the 
partition functions eqs. (|2.1|) and (|2.5j) . Such partition functions were already considered for one such 
pair in [H] when analyzing the quenched approximation and for n such pairs of degenerate mass in 
(Hj in the replica approach. The special feature of quarks and conjugate anti-quarks occurring in the 
partition function at the same time is that they may form a non vanishing meson density. This is, 
in turn, reflected in the existence of an effective chiral Lagrangian with isospin chemical potential as 
pointed out in linking it to QCD. 

Here, we will derive such an effective model in terms of a unitary group integral over the Goldstone 
manifold directly from the underlying random matrix model. We will allow for any, not necessarily 
equal number of quarks and conjugate anti-quarks with non degenerate mass. In the first model eq. 
(|2.1j) the presence of additional conjugate anti-quarks implies considerably more effort in computing 
the partition functions. In the procedure we will make use of the results from 6 . In contrast to that 
in the second model eq. 1)2.5)1 the resulting partition functions immediately follow from the theorem 
proved in 

The matrix model partition functions eq. 1)2.1)1 with m quarks of mass mj and n conjugate anti- 
quarks of mass 7ij in the sector of topological charge v is defined as 



Z 



(Nf=m+n,u) 



(n;{m f } m ,{n*} n ) = 

m I m/ljv i$ + /iljV \ n 

d$d& n det n det 



/=i 



n*l N 



\-i& + Lil ] N n* g l N+h 



(2.10) 



-N{qq) 2 TT3>& 



J 



Note that in contrast to the previous section we now allow for complex masses, as they may serve 
as source terms for the complex Dirac operator eigenvalues (see e.g. in 0). Using the standard 
fermionization technique we arrive at the following matrix model representation in terms of the complex 
matrix Q of size Nf x Nf 



dQdQi det[M + Q 



x det 



(M + Q f )(M + Q) - /x 2 (M + Q f )S 3 (M + Qt)-i Sg 



(2.11) 



N 



-N{qq) 2 TlQ^Q 



with the mass matrix M = diag(mi, . . . , m m , — n\, . . . , — n* ) and S3 =diag(l m , — 1„) being a general- 
ized Pauli matrix. The details of the derivation are given in appendixEI It is interesting here to work 
out the limit of real masses, M = Mt explicitly. To this end, we introduce a polar decomposition of 
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the matrix Q = U R with R = VrV^ being hermitian with positive eigenvalues r&, and U,V being 
unitary. Using eq. (|A.9|) and the Jacobian calculated in the appendix El we obtain 



zf'^foM) ~ e~ N ^ 2TlM2 f^f[dr k rl +l e- N ^l A Nf (r 2 ) 2 f dV J dU &et[tf] 



N 

Vr 2 V ] - u 2 VrVlWEzUVr- l VlY,z 



x e N{qq) 2 TlCM(UVfVUVfVW) ^ 

(2.12) 

This expression will be used in the next section. 

Let us now consider the second matrix model ^5] with pairs of complex conjugate quarks inserted. 
The corresponding partition function is defined as 

/N m n 

J] d 2 z 3 w( Zj ,z*) n rW}(z] + m 2 ) n n* g »(zf + O \ A n(z 2 )\ 2 
j=l f=l 9=1 

IN / N f n \ \ 

- (n n-/^+-/)n^f +o ) . 

\i=i \/=i 9 =1 / / 

where we again presented this object in a form of an expectation value. In contrast to eq. (|2.13|) we 
do not need to distinguish between signs of masses and chemical potential for the quarks and their 
conjugates. The reason is that the large-A?" result for eq. (|2.13|) turns out to be a quadratic function 
in all masses. Apart from that, in the identification between the nonhermiticity parameters r and /i 
appears quadratically. The same weight function eq. (|2.6|) for all eigenvalues is therefore equally valid 
for both signs of the chemical potential 3 . 

The evaluation of eq. ()2.13|) is again straightforward, due to the general theorem proved in |39| . To 
write the result in a compact form we need to introduce more notation. In addition to the orthogonal 
polynomials eq. (|2,8[) the so-called bare kernel made of these polynomials appears in the corresponding 
expressions: 

2V-1 

v 2 „,*2n, = pW^pCf),*^ 



k=0 

N-l 



_ 1 ^ T(u + l)k\ r 2k & )(Nz 2 \ {u) (Nu* 2 \ 

~ WHr)t Q WTkTTf [-*r) ■ (2 - 14) 

The kernel contains a sum over the orthonormalized polynomials pjf\z 2 ), 

p£\z 2 ) = h'h^Hz 2 ) , (2.15) 

with norms given by 

h t . J M^)Pf\:>)Pf\^) = /<"'M (l)' ^" 1 • (2.16) 



3 In this sense the model does not distinguish between isospin and baryon chemical potential. 
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Here we have introduced the functions 

/M(r) = Jd 2 zw(z,z*) = vrr^ + ^U-r 2 )^!^! (^=) , (2-17) 

where V^{x) stands for the Legendre function. The full kernel K^(z, u*) as it appears in the expression 
for the eigenvalue correlation functions is then obtained by multiplying the bare kernel with the weight 
functions, K n (z,u*) = [w(z, z*)w(u,u*)]z kn(z 2 ,u* 2 ). Following jHS] we immediately obtain 

B(mi,n*) 

where we may assume that n > m, without loss of generality (the case n < m follows from complex 
conjugation). This result is exact for finite- N. We note that the model eq. Q2.13JI (and eq. (|2.5jl l is 
always in the phase with broken chiral symmetry |18j . 

3 The large- N limit 

3.1 The weak nonhermiticity limit 

The limit of weak nonhermiticity ^B] is defined by taking simultaneously the large- iV limit N — > oo 
and the hermitian limit r — > 1 or fi — > in such a way that the following product is kept constant: 

lim iV(l-r 2 ) = a 2 (3.1) 

or, equivalently 

lim 2N(qq) 2 fi 2 = a 2 . (3.2) 

JV— »oo, fi— >0 

In this limit the macroscopic spectral density has support only on the real line and is given by a 
semicircle for both models eqs. (|2.1|) and (|2.5|) . in the limit r — > 1 and fi — > respectively .18 . In 
contrast to that the microscopic correlation functions differ from those on the real line. They still 
extend into the complex plane and depend explicitly on the parameter a or a. To identify them for 
the two models a relation between r and fi, or a and a, has to be imposed, which we will find by 
comparing the two partition functions. 

We will see that also the partition functions at weak nonhermiticity contain important information, 
and this has already been exploited for example in |]] when computing the microscopic density. It 
turns out that the partition functions with quarks alone differ from the partition functions at \x = 
only by an overall prefactor exp[— Nja/2]. When adding conjugate anti-quarks the situations however 
changes and the partition functions differ from their fi = values nontrivially 4 . 

The microscopic rescaling of the quark masses (and complex eigenvalues) at weak nonhermiticity 
is defined as 

Cf = 2N(qq)m f ,/ = 1, . . . ,m 

C g = 2N(qq)n* g ,g = l,...,n (3.3) 

4 Obviously, for only conjugate anti-quarks alone we are back to the situation of only quarks. 



n 



f=i m f 



n 



n 

9=1 



A m (m 2 )A n (n* 2 ) 




det [B(mi,n*)} , (2.18) 
i,j=l,...,n 



P 



-m 2 , -n* 2 ) 

*2\ 



N 



1, . . . , m 
m + 1, . . . , n 



(2.19) 
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where 2N(qq) ~ p(0) is the macroscopic density at the origin of the model eq. (|2.1|) . after taking the 
limit [A — » 0. 

We begin with the partition functions containing only quarks from subsection 12. 11 First of all due 
to the rescaling of the masses eq. 1)3.3)1 the constant prefactor e-^gg^Tr Af 2 reduces to unity. Taking 
the weak limit eq. ()3.2j) we can replace the factor 



(H ~ M ) 



2\N 



J2N 
r k 



exp [-Nil 2 



(3.4) 



inside the integral eq. 1)2.4)) ° ver the radial coordinates. After rescaling the arguments of the Bessel 
functions through eq. ()3.3|) a saddle point evaluation in the variables r& leads to the value 



^Isp = lNf(qq) 1 



(3.5) 



Since the determinant in eq. (|2.4j) is degenerate at the saddle point we have to take into account the 
Gaussian fluctuations. Performing the calculation, we get 



zf'^MC/Dlweak ~ ™p[-N f {qq) 2 N^ 



detjfc. 



i=i, 



A^v f (C 2 ) 



(3.6) 



With li j \x) we denote the j— th derivative of the modified Bessel functions. This final expression for 
the partition function generalizes previous results which were known only in the sector of topological 
charge v = with Nf = 1 ^2] an d Nf = 2,3 flavors 15 . Eq. (|3,6|) is exactly the same as for 
\i = [HI EH]) apart from the exponential prefactor exp[— Nfa 2 /2\. This implies in particular that 
the partition functions obey the same consistency conditions |44| as those for /i = 0. More generally 
speaking, they belong to the same Toda lattice hierarchy as it was already exploited in 



We turn to the second partition function eq. 1)2.9)1 . In order to read off the proper rescaling defined 
in eq. 1)3.3)1 we first take the Hermitian limit r — > 1 of the weight eq. (|2.fij) in order to determine 
the corresponding variance. This procedure is based on the macroscopic spectral density at weak 
nonhermiticity being given by taking the Hermitian limit 5 : 



lim w(z, z*) 

T— >1 



lim (x 2 + y 2 Y + z exp 



T-fl 



„2v+l 



N 
1 + t" 



N 



-x 



1 - T 



y 



7rcr 
2iV2 



5 (y) exp 



(3.7) 



where we have inserted z = x + iy. Therefore we have to use the value (qq) = l/\/2 in eq. (|3.3|) . 
leading to the correspondence Q = V^Nmj and similar for the conjugate. Such a rescaling which is 
different for the two models is precisely the mapping between the two sets of different mass parameters 
mentioned after eq. (|2.6[) . 

We could have also introduced a variable (qq) in the model eq. (|2.9j) by rescaling the eigenvalues. 
However, since the parameter drops out after microscopic rescaling we kept (qq) 2 = 1/2 for simplicity 
here. In order to perform the weak nonhermiticity limit in eq. 1)2. 9|) we first extract the powers 
of r used in the definition of the polynomials 1)2.8)1 from the corresponding determinant. Using the 
definition r 2 = 1 — a 2 /N we obtain the prefactor 



lim 

N^OO T- 



T 



N f N+N f (N f -l)/2 



>1 



exp 



-Nfa 



(3.8) 



The corresponding weak nonhermiticity limit for the weight function is given in section Hj. 31 eq. 13.591 . but we do 
not need it here. 
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The remaining determinant of Laguerre polynomials can be evaluated in the same way as for [i = 
using eq. (|3.3|) 

-Nm 2 \ 

UO , (3-9) 



lim m u L^ 



AT— KX),T— >1 



J N 



2t 



and differentiating due to degeneracy. Ignoring constant factors and powers of N we thus obtain 

l detfcj = i v .. j Ar / 



Z 



(N f ,u) 



11 



«;{C/})I 



weak 



exp 



--N f a 2 
2 ; 



(3.10) 



A^(C 2 ) 

We immediately see that the two partition functions ()3.6|) and H3.10|) agree upon identifying 

(1-r 2 ) = 2<gg)V (3.11) 

or, equivalently a = a. Such an identification can be conveniently interpreted in terms of equating 
the two macroscopic densities in the complex plane for small \x ^Hj. In principle, partition functions 
can always be multiplied by an overall constant, making them agree. The nontrivial statement is that 
in both cases the rescaled nonhermiticity parameter factors out. 

We thus proved the equivalence of the two partition functions is the limit of weak nonhermiticity, 
for an arbitrary number of flavors Nf of nondegenerate masses mf for any given topological sector v. 
This equivalence was claimed in |15| and verified there for degenerate masses and sector v = 0, up to 
and including Nf = 3 flavors. 

Next we investigate the weak nonhermiticity limit of the partition function with quarks and con- 
jugate anti-quarks, eq. (|2.11|) . In this case, all powers higher than linear in M and fi 2 in eq. (|2.11j) 
are dropped. After expanding also the determinant to the power TV" in powers of M, we have 



( M ;M) 



dQdQi det[QTe 



(3.12) 



x exp 



iVTrln(QtQ) + _/VTr (m^ 1 + Qt-1) _ f^Q-^Q^ 1 ^ ) + 0(1/N) 



By using again the polar decomposition Q = U R and the respective Jacobian, we obtain 



J JI dr * rl +l <r N ^ 2r l +m * r l A Nf (r 2 ) 2 jdVjdU det 
exp [iVTr (^{Vr^V^tf + UVr^V^)) - /x 2 iVTr (vf-^U^UVf^V^ 



(3.13) 



In this form the integral is amenable to a saddle point approximation at large- N in the variables r^, 
which will lead to the chiral Lagrangian picture In the previous case, with only quarks included 
the /^-dependence factorizes out as we discussed above. Now the dependence will be less trivial, making 
it a real check for the equivalence of the two models. Due to the rescaling eqs. 1)3 .3(1 and (|3.2jl the 
exponents inside the unitary integral are of the order of unity and the the saddle point value is taken 
at 

r| sp = lNfiqqy 1 , (3.14) 
as previously in eq. ()3.5|) . Since the unitary integral is non vanishing at r| S p ~ ljv>, we have 



sp 



dl7 det [E^]" exp N(qq)Tr(M(rf + U)) - fi 2 N{qq) 2 Tr(U^ 3 UE d 



(3.15) 
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This is precisely (a zero-dimensional version of) the chiral effective theory for isospin chemical potential 
as derived in with replacements 2N < > V and fJ- 2 {qq) 2 <- 



— > H F*. The manifestation of the 
zero-dimensional nature is that the parameter F n , the pion decay constant entering in the chiral 
Lagrangian is not contained in the present matrix model partition function, as it comes from 
the non-zero momentum modes of the Goldstone bosons. If we wish to write the [i— dependent term 
as iTr[Z7', S3] [U, S3] as it appears when deriving it from a gauge principle, we have to multiply a 
factor exp [N (qq) 2 /j? N to the partition function. By doing so also in the previous section it would 
remove the trivial //—dependence, that leads to an unphysical negative quark number density in the 
broken phase (see e.g. in ^2] f° r Nf = 1). Finally we wish to mention that the usual mass term in the 
unquenched, effective QCD Lagrangian is Tr(M^U^ + MU) instead. There complex masses are needed 
to locate the dependence of the 9— angle in the action for example, leaving the partition function real. 

We are now ready to evaluate the unitary integral in eq. (j3.15j) exactly, following closely There, 
the following parameterization of the unitary matrix U G U (Nf = m + n) has been suggested 6 



U 



A 



vi 

v 2 



m 
u 2 



A 



/ A Vim ~ A 2 

Vim ~ A 2 "A 

V 



v[ 
vt 






-ln- 



(3.16) 
(3.17) 



Here, we denote A = diag(Ai, . . . , X m ) with Xk G [0, 1] for k = 1, . . . ,<m. The unitary submatrices are 
Ui,v% G U(m), U2 G U(n) and v 2 G U(n) = U(n)/(U(l) m x U(n — m)), where n>m has been chosen 
without loss of generality. The Jacobian for this transformation was computed in [5], 



J{{^k}) = n(Ag - A?) 2 Yl2X k (l-X 2 ) 

k>l j=l 



(3.18) 



For the trace containing S3 it follows 

Tr(C/ts 3 c/S3) 



n 



3m + 4^Al . 



(3.19) 



k=l 



With these steps taken we can insert the above parameterization in eq. (|3.15|) . arriving at the following 
factorized group integrals 



rl m 

Zl Nf (M; {™>f} m , {n* g } n ) | weak ~ exp [-(n - 3m)N(qq) 2 [i 2 ] / TT d\ k J({A fc })e 



-A(qq)*Nl?\l 



I du\ I dv\ det^i]" exp 

IU{m) Ju(m) 

x / dv,2 dv 2 det[u2] l/ exp 

JU(n) Ju(n) 



N(qq)Tr(u\v \rh v\X + v\rh v\U\X) 



N (qq)Ti(u\v2n*v 2 X- + v\n*v 2 u 2 X 



(3.20) 



By splitting into blocks we introduced the following obvious notation, denoting rh = diag(mi . . . , m m ), 
h* = diag(— n\ . . . , — n*) ; and by A_ = diag(— Ai, . . . , — A m , —1, . . . , —1) a matrix of size n x n. Per- 
forming the additional transformations U{ — > v]ui and renaming m — > u] for i = 1,2, we can bring 



3 In addition to we have performed an extra rotation u% — > Vimvj for i = 1, 2, interchanging the first 2 matrices. 
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both double unitary integrals to the form 



du / dv det [u v] v exp 
u Ju 



— Ti{u^d v^b + v a ub) 



——teMoibj)], (3.21) 



with U = U(l) for / = m or n, respectively 7 . Here we have used the result 45_ which gives this integral 
also in the case when the matrices a and b are not hermitian. In our special case the af and bf are 
the complex eigenvalues of the diagonal matrices d 2 and b respectively. In the second double unitary 
group integral we still have to take limits due to the (n — m)-fold degeneracy in the matrix A_. Both 
integrals together then cancel the Jacobian J({A^}) up to the factor HfeLo^Afc, leading to 



^WC/WWkveak 



exp[— (n — 3m)(qq) 2 Nfi z 



*2\ 



x det [I„(A fc Cl)] det 

k,l=l,...,m 



A m (c 2 )A n (e 

/ / v (ejAi) ••• Iv&Xm) /„(£) 



I m 



fee 



-4{qq) 2 N^X 2 k 



k=X 



t*n— ro-1 r(n-m-l) 
*n-m-l r(n-m-l) 



(3.22) 



The above equation can be further simplified. First we multiply the two determinants of Bessel 
functions, det (5) det(yl) = det(AB). Because of the different size of the two matrices we have to add 
a block of unity l n -m to the first matrix Iv(\kCl) °f si ze mX in. Each entry in the first m columns 
of the resulting product matrix AB contains a sum of m terms, Y^j=x -^(^&^')^(^iCi)- Due to the 
symmetry of the integrand under permutations of the variables A& and due to the invariance properties 
of determinants the sums can be reduced to single terms, I v (^Xk)I v (XkCl) i an d we only sketch the 
procedure briefly. We begin with the first column. As determinants differing only by one column can 
be added, we expand the determinant as a sum of m terms, which only differ from each other by the 
labeling of A& in the first column. Due to the invariance with respect to permuting A/% all terms can be 
written as m times the same determinant, with label Ai in the first column. Next we can eliminate all 
the terms with label Ai in the remaining columns, by successively subtracting the first column times 
an appropriate factor. Next we process the second column in the same way, keeping only the label A2 
there, and so forth. 

We finally arrive at the following expression 



Z 



(N f=m+n,v) 



(MC/}m,{C}«)lweak 



exp [— (n — 3m)(qq) 2 N/j, z 



A m (c 2 )A n (e 2 ) 

/ HdX 1 X 1 e-* N ^)'^I v ^X 1 )I v (X 1 Ci) ••• J 1 dA 1 A 1 e- 4Jv ^) 2 ^ A ?J 1/ (CA 1 )I,(A 1 Ci) \ 



det 



Jo 



dX X e~ ml - qq i^ x ^I (tfX )I ( A C ) 



Jo dX m X r 



n—m- 



1 Jn-m-X) 



)Iu (AmCr? 



t*n-m-X An-m-X) 



(3.23) 



7 In order to apply the result |45| we have to promote V2 to the full unitary group U(n). We first multiply by the 
additional integrations (f u ^ dwi det(wi) v ) m f [f ^ n _ m ^ dwo det(wo) L ' and shift U2 — >diag(n;i , . . . , w m , wq)u2. Due to the 

cyclicity of the trace and the (n — m)— fold degeneracy of A_ which makes it commute with diag(^i, . . . ,w m ,wo) we 
obtain a matrix U2diag(wi, . . . ,w m ,wo) that parameterizes the full group U(n). 
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suppressing the symmetry factor ml as well as taking the transpose of the matrix. At the last step 
the m integrations over the have been taken inside the rows of the determinant. Eq. ()3.23j) 
generalizes the results obtained in [Hj for an equal number m = nof quarks and conjugate anti-quarks 
of degenerate, complex conjugate mass rh = zl n and n* = z*l n each. In jHj the partition function is 
given as well for n > m, for real degenerate masses rh = xl m and n* = yl n - 

We can now compare to the second model by taking the weak nonhermiticity limit of eq. (|2,18|) , 
To achieve this we first multiply all powers into the first m rows and all powers n* u into all n 
columns of the determinant B(rrii,n*). The large-iV limit of the norms eq. ()2.16j) and the bare 
kernel eq. (|2.14|) have already been taken in [ISj. Ignoring all factors of N the product of norms will 
lead to 

N+m-l 

lim TT h k ~ a m . (3.24) 

k=N 

The weak nonhermiticity limit of the bare kernel is given by ^H] 



lim m^n* v K/vh 



i(-m 2 



-n 



*2 



a 



d\\e~ a x * I u (\( f )I u (\C 



(3.25) 



where we have continued to negative arguments of the kernel. We have used again the value (qq) = 
l/v2 in eq. (|3.3|) corresponding to our model. Since all matrix elements of kernels have different 
arguments there is no degeneracy in this part of the determinant. Multiplying with eq. (|3.24|) cancels 
all inverse powers of a from the kernel inside the determinant. The part of the determinant containing 
only polynomials P^ u \—n* 2 ) can be dealt with as previously, using eqs. (|3.8|) and (|3,9|) . First we take 
out all powers of r, leading to 



lim T {n-m)(2N+n+m-l)\ ^ exp 

JV^OOT^l 



in 



m)a 



(3.26) 



The rows of polynomials become degenerate after taking the limit eq. ()3.9j) . leading again to differen- 
tiations as in eq. H3.1U|) . Performing the manipulations, we arrive at the following result 



r (N f=m+n,v) 



(a;{C/}m,{£g}n)| wea k 



exp [— |(n 



x det 



* n-m-l r (n—m—l) 



A m (c 2 )A n (e 2 ) 

•• fidXXe-^UXCMXZ*) \ 



ttdXXi 



i v {xu)i v {xti 



1 r (n-m-l) 



(Cn) J 



(3.27) 



The equivalence of the two partition functions can be established as follows. If we multiply eq. (|3.27|) 
by the overall constant expf+^A^a 2 ], with Nf = n + m, and then identify 



2 ..2 



(1-r 2 ) = 4(gg)> 



(3.28) 



or equivalently a = 2a the two partition functions ()3.23j) and ()3.27|) agree. We conclude that both 
models can be mapped onto each other as a function of rescaled masses and rescaled nonhermiticity 
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parameter, for an arbitrary number of quarks and conjugate anti-quarks of different, nondegenerate 
masses. In the present case the dependence on a is much less trivial as it enters all the integrals 
inside the determinant. We are not able to provide a simple explanation for the fact that the two 
identifications eqs. (|3.11|) and (|3.28f) differ only by a factor of two. However, our main point here is 
that there exists a mapping of parameters making the two partition functions the same. We will find 
yet another mapping of parameters in the limit of strong nonhermiticity below. 



3.2 The strong nonhermiticity limit 

The limit of strong nonhermiticity is defined by keeping fj? or r 6 [0, 1] fixed, independent of N when 
performing the large-N limit. Furthermore, the quark masses (and the complex Dirac eigenvalues) 
have to be rescaled with a different power in N |18l I15j 8 



, m 



Cf = V2N(qq)m f J 

C g = V2N(qq)n* g ,g = l,...,n. (3.29) 

In this limit complex eigenvalues of the Dirac operator fill in a truly two dimensional domain in the 
complex plane and the corresponding spectral density is nonvanishing there. With these definitions 
at hand we can again make a saddle point approximation for the partition functions in question. 

We start with the first model eq. (|2.4jl which contains only quarks. Due to the rescaling of the 
masses the arguments of the Bessel functions now become large. But even when replacing them 
with their asymptotic value, I u (x) ~ exp[x]/y/x, they will not contribute to the saddle point, as the 
argument is of order 0{yN) and is small compared to terms 0{N) in the exponent for radial variables. 
The saddle point is thus given by conditions 

2r k 



and we choose the positive value 



0, k = l, 



f 



(3.30) 



(3.31) 



Hsp = J- AT, 

as the relevant solution. For a more detailed analysis, in particular, on the connection of the saddle 
point solution at zero to the phase transition we refer to section |I] below. Taking the usual limit of a 
degenerate matrix at r|sp ~ lzv> we arrive at 



4 JV/ ' iy) (^;{C/})l s trong ~ exp[-NN f (l + (qq) V)] exp 



N f 

■lEc 



k=l 



det 



k,l=l,...,Nf 



>fc-l r (fc-l) 



^QV 1 + (w> V 



Atv^C 2 ) 



(3.32) 



Now we treat the second model by looking at eq. (|2.9|) at finite- N. The strong nonhermiticity limit 
is most easily taken here by using an integral representation of the Laguerre polynomials, following 
IE] for N f = 1, 



ds e 



-Ns g N+Z 



2j„(2VsNx) 



X 2 - 

e e 



■N 



J u (2VNx) 



(3.33) 



For the constant proportionality factor we have kept (qq) as in the weak limit, although the macroscopic spectral 
density will in general no longer be constant but rather depend on /i. For a fj,— dependent rescaling we refer to the 
discussion after eq. (13.361 and after eq. 13,4711 . 
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At the second step we have made a saddle point approximation, taking into account that x 2 = 
—Nm 2 /(2T) is fixed and finite. Using the rescaling of the masses eq. ()3.29j) with the chosen value 
(qq) = 1/V2, corresponding to yNm = (, and taking care of the degeneracy of the determinant we 
obtain 



^f /,l V;{C/})lstrong ~ exp[-NN f ]r 

detfc ; = i w 



N f (N+i(N f -l)+±v) 



27^ Cfc2 

k=l 



(f 1 ^ (V2NC1T-* 



Using now the relation eq. 



Aiv f (C 2 ) 



r = y/1 - 2<gg) V 



(3.34) 



(3.35) 



we see that the //—dependence does not match the previous case in general, unless we expand for small 
ji 2 . To the leading order we obtain 



^'^(MC/Dlstrong ~ exp[-JViV>(l + (g g )V 2 )+O0i 4 )]exp 



det fc 



l=l,...,N f 



C k-lr 



JV f 

^(i + <w>V + o(a* 4 ))Ec2 

2Aud +i(??)V + 0(M 4 ))) 



fc=l 



(3.36) 



The two expressions eqs. (|3.32|) and (|3.36|) match to the leading order terms in the determinant and 
in the first factor. However, the coupling between masses Cfc and [i introduced in the second model 
in the second exponential prefactor does not have an analogue in the first model. Here the matching 
only holds to terms of zeroth order. 

It is worth mentioning that even apart from the //—dependent exponential suppression factor 
exp[— NNf(qq) 2 n 2 ] a proper large-iV limit of the partition functions in terms of rescaled masses 
does not exist. The Bessel functions still depend on \/~N in the argument. In the asymptotic limit 
however, the exponential suppression wins. Our microscopic rescaling of the masses eq. (j3.29|) cannot 
be modified to achieve an iV-independent result. Such a rescaling is, in fact, dictated by corresponding 
rescaling of the complex eigenvalues ^H] necessary to obtain a smooth limiting eigenvalue correlation 
functions. This type of rescaling in N at the regime of strong nonhermiticity was also found earlier 
in models without chiral symmetry |38j . 

One also may wonder if it is possible to match the two partition functions eqs. (|3.32|) and (|3.34j) 
beyond an expansion in /t, by introducing a different, /x- and r-dependent rescaling of the masses 
in eq. (|3.29|) respectively. However, a little thought shows that due to essentially different r- and 
//-dependence of the masses of the two models such a procedure is impossible. 

In order to take the large- N limit at strong nonhermiticity of the first model with quarks and 
conjugate anti-quarks we go back to the expressions for finite- N of the partition function, eqs. ()2.11j) 
At the regime of weak nonhermiticity it was crucial to expand the determinant det 



and (I2TT2 

to be able to perform exactly the unitary group integrals arising from the parameterization of Q. 
However, at the regime of strong nonhermiticity such an expansion is no longer possible with respect 
to the parameter /i since the latter no longer scales with N. For the sake of simplicity we will restrict 
ourselves to real quark and conjugate anti-quarks masses, M = M\ eq. (|2.12|) . In general the radial 
and unitary degrees of freedom no longer decouple. However, if we assume that the permutation 



N 
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symmetry makes all radial coordinates to take the same saddle point value as it was the case 
previously, see eq. (|3.5j) . that means 

^Isp = Tspl^ , (3.37) 

with the value r S p to be determined, the integral eq. ()2.12[) simplifies considerably. The unitary 
matrix V drops out from the integrand and we obtain 



zf^V;M) ~ e" 



N {qq) ^TM^-N(m?N f rl v f ^^t] 



x exp 



N(qq) 2 r sp TrM(U + tf) 



det 



^sp 1 ^/ 



2 C/ t E 3 l7Sf 



N 



(3.38) 



This expression can be regarded as the effective partition function at the regime of strong nonher- 
miticity, and serves as generalization of the effective Lagrangian from chiral perturbation theory 
eq. (|3.15|) which was valid at weak nonhermiticity. The two expressions only agree after expanding 
to the first order in [i 2 (for real masses). Note that while in the chiral Lagrangian picture higher 
powers terms (^U^^UY,^ can be excluded from power counting in U, they cannot be excluded in 
the matrix model. We will see that in order to describe a phase transition to the symmetric phase it 
will be important to keep all powers in fj,, see section below, also cf.|12|. 

In the following we also restrict ourselves to an equal number of quarks and their conjugate 
partners, n = m. The reason is that in the general case n > m one naturally expects a result to 
be of the mixed form (compare eq. Q2.18JI or eq. (|3.23|) . containing both limiting "kernels" and 
"polynomials" 9 . However, we have just seen that the quantities corresponding to the polynomials do 
not have a proper large- N limit. Since we are looking for partition functions that do possess such a 
limit at strong nonhermiticity (as we will find for the second model below) we choose n = m to ensure 
that the result contains only kernel terms. 

In eq. 1)3.38(1 we can again employ the parameterization eqs. 1)3.16)) and (|3.17|) . where it is instruc- 
tive to first look at the determinant alone. We obtain 



det 



sp 



N 



det 



det 



r spl" 



/uV(2A 2 



2/j, 2 v 1 \y/l n - A 2 )? 



N 



r 2 1 

' sp A r 



r 4 1 



2/zWln - A>1 

2r| p /x 2 (2A 2 -l„) + /i 4 l 
" ' 4r 2 p ^ 2 , 



/j, 2 v 2 (2X 2 - l n )v. 



N 



N 



n i 



i=i 



(r 2 S p+^ 2 1 



(3.39) 



where all angular-variable dependence has dropped out. The integrals over the unitary subgroups ui t 2 
and 1)1,2 of the mass dependent exponential ex.p[N(qq} 2 r S pTrM(U + U')] can be performed as in the 
weak nonhermiticity limit eq. (|3.2()j) . using the integral eq. (|3,21|) , In fact we can literally repeat all 
the following simplifying steps there after eq. ()3.22j) . leading to 



2f^V;{C/}n,{£ fl }r 



A„(C 2 )A n (e 2 ) 



-iy 
e 2 ^ 



k=1 «l+e k ) r Nf(-+l) p -N(m) 2 N f rl p (r 2 2) 2nJV 



sp 



sp 



x det 

i,j=l,...,n 



d\X 1 



4r 2 p ^x 2 



N 



(r 2 



sp 



2\2' 



2N{qq)r sp XCi)lu [V2N(qq)r sp XZ j 



(3.40) 



Here we have also inserted the microscopic scaling of the masses eq. (|3.29|) . We still have to determine 
the saddle point value r S p and take the large- iV limit. In the latter the integral can be computed by 

9 For the interpretation of the partition function in terms of kernels and polynomials we refer to section EQ1 
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making the change of variables s = 4rgp/i NX j (r$p + fi ) and using that limjv^oo(l — s/N) 



exp[— s]. We also employ the following integral, 

J dse~ s J v (V2sC) J» (V^C) = exp 



urn 



(3.41) 



after analytically continuing in the masses 10 . Taking out common factors of the determinant we arrive 
at 



^V{C/}n,{U«) 



1 



A n (C 2 )A n (£ 2 



e 2 



^ n k=1 (C 2 k +e k ) e - N ^ 2N f r2 S p ^ + ^2nN 



e 2>™ ^'SP^ I ^k=l ~^T- 



det 

i,j=l,...,n 



4/i 2 



The value for rsp can finally be read off as 
which allows to simplify the final result down to 



-2 



i 



A n (C 2 )A„(£ 2 ) (4(gg) 2 /i 



-2nN{{qq) 2 ^-l) 



(3.42) 



(3.43) 



(3.44) 



x exp 



J(i-4(^)V)El^ f) 



\ 4(<?<?)V J U 



det 

=l,...,n 



7 r ^ 



Our conclusion is therefore that the large-iV limit of the partition function (as a functions of the 
masses) with an equal number of quarks and conjugate anti-quarks is well-defined, in contrast to that 
for quarks alone, eq. (|3.32j) . The /x-dependence can be almost entirely absorbed by redefining the 
rescaling of the masses eq. (|3.29|) to Q — > Q /(2(qq){i) and £ 9 — » (, g /(2(qq)fi), giving: 



^'""Vi {C/}», {&}») 



-2nN{{qq) 2 ^ 2 -l) 



A n (C 2 )A n (£2) (4<gg)V) n 



e i(i-^W>V)E2 = i«2+€» det 



i,j=l,...,n 



The strong nonhermiticity limit of the second model with quarks and conjugate quarks is again 
performed easily, without being restricted to real masses. Looking at the expression eq. (|2.18|) we 
need the limiting normalization factors, kernel and polynomials. The latter have been analyzed by us 
already, and we start with the normalization factors. Omitting all constant factors we obtain 



JV+771-1 

n 

k=N 



(f (u) (r)Y 



(3.45) 



which contains all r— dependence. Before taking the strong nonhermiticity limit for the kernel (and 
polynomials) we multiply all the prefactors and n* u into the determinant det[B(m,i, n*)] in eq. 
(|2.18|) . as in the weak limit before. The asymptotic kernel has been already evaluated in |18j . reading 



v * v / 2 

m n KN[—m 



-n 



,*2\ 



(1 



T 



/M(r) 



exp 



(C+r 2 ) 



2(1 -r 2 ) 



i 



r 



(3.46) 



10 The integral in eq. 13.4011 is also convergent after change of variables and taking the large-TV limit as it stands, 
without continuing to imaginary masses. 
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Multiplying all factors of f^ v \r) from eq. (|3.45j) into the first m rows of the determinant cancels the 
(r)-dependence of the kernel. Using eq. (|3,33|) for the asymptotic polynomials results in the full 
expression for the partition function reading 



7 (N f =m+ro,f ) 
J II 



x det 



( r ;{C/}m,{£g}n)lstrong " 



^—N(n—m) 



_(n-m)(N+ \ (n+3m~l+2u)) 



A m (C 2 )A n (e 2 ) (1-t2)> 



e; 2 



i>-m—lj(n—m—l). 



,* 2 
^71 



2t -^*n-"»-l/( rl - m - 1 )( v /2iV^* r -2) / 



(3.47) 



In order to compare with eq. (|3.45|) we look at the particular case of equal number of quarks and 
conjugate anti-quarks, n = m, 



^II f ^( T i {C/}n> {Cg}n.)lstrong 



A n (C 2 )A n (e 2 ) (l-r 2 )« 



exp 



x det 

i,j=l,...,n 



1-r 2 



(3.48) 



This expression has a finite limit at large TV as a function of r and the masses. We note that 
here the masses always appear in the same r— dependent combination l/Vl — t 2 . If we define a 
t— dependent microscopic rescaling my/2N/(l — r 2 ) = £ instead of eq. (|3,29|) that square-root factor 
can be absorbed: 



^■11 f ( r ; {Cf}n, {Cg}ra)lstrong 



t (c 2 )A n (e 2 ) (i - r 2 )" ij=t,n [lu Ui Cj jJ • 

(3.49) 



We note that p(0) = 2^(i3p) * s ^ e cons tant macroscopic density corresponding to uniform filling of 
an ellipse. The scaling defined above is very reminiscent to the scaling on the real line used to get read 
of the mean level spacing. The two partition functions eq. Q3.45JI and (|3.49|) can now be identified 
upon the mapping 

t = 1 - 4(qq) 2 fi 2 . (3.50) 

It is yet different from the two previous mappings, the main point of our analysis being the very 
existence of a correspondence between the two functions in terms of their arguments. 



3.3 Universality 

After having demonstrated several cases of agreement between the two partition function with a given 
flavors content let us use those results for discussing the issue of universality of the random matrix 
theory in the complex plane. For matrix models with eigenvalues on the real line universality is a 
well established principle, both on the base of heuristic, physical arguments |21l 1221 136j as well on 
firm mathematical grounds (see |461 133j and references therein; we mainly discuss works relevant for 
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QCD applications). On the other hand, for complex eigenvalues the subject has so far attracted only 
little attention, apart from the discussion in [23]. As we already have mentioned several times, it 
is important to distinguish between the limits of weak and strong nonhermiticity ^Hj. As is well 
known, the eigenvalue correlations in the weak nonhermiticity limit interpolate between those typical 
for real eigenvalues and those known for complex eigenvalues. More precisely, sending the rescaled 
nonhermiticity parameter a from eq. (|3.2|) to or oo allows to recover correlation functions typical 
for eigenvalues of Hermitian or, respectively, complex strongly nonhermitian (Ginibre-like) matrices. 
This was explicitly checked in several cases, as for example in ^H] for the unitary (UE) and in ^B] for 
the chiral unitary ensemble (chUE). 

In that sense the weak nonhermiticity limit is closely related to real, universal correlations, and it 
is natural to expect it should maintain at least some universality features. In [23] a certain class of 
deformations of measures of the complex UE were studied with the conclusion, that the asymptotic 
polynomials, kernels and correlation functions remain universal. 

It therefore came as a surprise that for the complex extension of the chUE two slightly different 
results were obtained for the microscopic spectral density. Indeed, the second model eq. Q2.5|) with 
complex eigenvalues was solved in the paper ^H] by analytical continuation of the Laguerre polynomi- 
als, and a comparison to quenched QCD lattice data with chemical potential confirmed the obtained 
results |19j . On the other hand, the original model eq. (|2.1[) based on a random matrix representation 
as a starting point was solved recently for the microscopic spectral density in the weak limit, using an 
exact replica approach, see UJ. The obtained density, although having a very similar structure, was 
slightly different, and the difference was confirmed by simulating numerically the underlying random 
matrix model, eq. (|2.1[) . We remark however, that for the range of parameters \x used in the analysis 
of the lattice data Jlj the difference was too small to be detected. In view of the agreement between 
partition functions of the two models found in the preceding sections and in ^Sj we will try to shed 
some more light on the origin of the discrepancy, and thus on the issue of universality. 

In the random matrix theory several different objects can be tested for the property of universality. 
Let us start with recalling some basic facts for models that can be solved with the technique of 
orthogonal polynomials, such as eq. IJ2.5|) or the chUE. The solution for the eigenvalue correlation 
functions can be written as 

k 

p k {z!, . . . , z k ) = det [K N (zi,z*)} = Y\w(zi,z*) det [K N (zi,z*)} 

i,j=l,...,k J i,j=l,...,k 



Y\w(zi,z*) det 



1=1 

N-l 



fc=0 



(3.51) 



where we use the notation introduced in eqs. (|2.14|) and (|2.15jl . 

It is evident, that three different objects can be analyzed from the point of view of asymptotic 
universality in the large-iV limit: the polynomials Pj^\z), the bare kernel Kiy(zi,Zj) and the weight 
function w(z,z*) itself. All these objects have a direct relation to partition functions as we will see 
below. For the universality of the spectral correlations eq. (|3.51l) we need that both the weight and 
the bare kernel are universal, as the full kernel is given by Kjsr(zi, z*) = [uu(zi, z*)w(zj, z*)] 2 kn(z{, z*). 
We also note that for real eigenvalues the universality of the bare kernel follows directly from that of 
the polynomials 11 due to the Christoffel-Darboux formula 

K N (x,y) = hZ 1 , . (3.52) 

x — y 



L We note however, that the universality of the norms hk is a separate issue, see e.g. in |21|. 
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For orthogonal polynomials in the complex plane this relation does not hold and thus the universality 
of polynomials does not necessarily imply the universality of the bare kernel, and vise versa. 

There exists an alternative method for computing spectral correlation functions from working with 
characteristic polynomials. An advantage of the method is that it retains its validity also in a general 
case when orthogonal polynomials are unavailable, for example, when the partition function cannot 
be represented in terms of eigenvalues. In the simplest case the resolvent G(z) can be generated from 
differentiating a single ratio of characteristic polynomials, 



G(z) = (TV 



1 



z-H 



N 



detjz - H) 
det(tt - H) 



N 



(3.53) 



Here H = denotes a hermitian random matrix of size N x N averaged over a Gaussian or more 
general weight function. Knowing the resolvent, the spectral density follows by taking the discontinuity 
along the support or the antiholomorphic derivative, for real or complex eigenvalues, respectively. For 
getting access to higher order correlation functions more ratios of characteristic polynomials can be 
used as source terms. 

Apart from being a generator for spectral correlation functions characteristic polynomials can be 
regarded as interesting objects in their own right, and it is therefore natural to ask about their uni- 
versality. In fact, it is possible to relate characteristic polynomials to orthogonal polynomials in such 
models where the latter are known. First of all, the expectation value of a single characteristic poly- 
nomial directly gives the orthogonal polynomial for the corresponding model in monic normalization: 



(det(z - H)) 



N 



Pn(z) 



(3.54) 



Second, the expectation value of the product of two characteristic polynomials yields directly the bare 
kernel (see e.g. in j47j ) 



(det(z - H) det(u - H)) 



N 



h N k n+ i(z,u) 



(3.55) 



Third, for real eigenvalues the expressions involving the inverse of a characteristic polynomial are 
given in terms of the Cauchy transform of the orthogonal polynomials, as was recently discovered in 
|34"j . and further developed in [HH1 E""3 E21 IHEj ■ For example: 



det(z - H) 



N 



_L f dx^-P N {x) = Hz) , 

Z7TI / Z — X 



(3.56) 



where to ensure that the object is well-defined, the poles have to be suitably avoided by giving an 
imaginary part to z, 9m(z) ^ 0. Various formulas expressing arbitrary products [271 I2SI EE1 EE9 
and arbitrary ratios |341 1331 13*5] of characteristic polynomials in terms of orthogonal polynomials, their 
Cauchy transforms and bare kernels containing one or both of these have been proved recently. 

A natural question which immediately arises is whether the universality also holds for arbitrary 
ratios of characteristic polynomials, especially in view of their relation to orthogonal polynomials 
and bare kernels just mentioned. We would like to point out that because of the presence of the 
Cauchy transform the universality of such objects is in general not a simple consequence of the known 
universality of the kernels and polynomials. For the unitary ensembles with real eigenvalues this 
question has been completely answered. For the standard UE the universality was rigorously proved 
in |33j . Arbitrary ratios of characteristic polynomials for the chUE have been computed in |301 1321 l3*T] 
and proved to be universal in various regimes in |36| 13 7j . Let us mention that all what we have 
said immediately applies to the matrix model partition functions of QCD, as the insertion of massive 
flavors is nothing else than the insertion of Nf characteristic polynomials. In fact, from the point of 
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view of massive partition functions the universality of arbitrary products has been proved previously 
in H2] for the UE and chUE. 

For characteristic polynomials of matrices with complex eigenvalues much less is known. First 
of all we have to distinguish between the characteristic polynomial and its complex conjugate, just 
as we have distinguished between quarks and conjugate anti-quarks in our previous considerations. 
For arbitrary products of characteristic polynomials and a different number of complex conjugate 
characteristic polynomials an expression in terms of polynomials and bare kernels has been derived 
in |39| (see eq. (J2.18|) which we have already used). From what has been said in this section it is 
clear, that using products alone we will not be able to deduce the universality complex eigenvalue 
correlations. 

We can thus conclude the following. We have studied two different models, eqs. (j2.1|) and (|2.5|) . 
which are apriori not the same, in particular as the former does not admit a simple eigenvalue repre- 
sentation. We have found that in the large- ./V limit at weak nonhermiticity the two partition functions 
agree for an arbitrary and different number of quarks and conjugate anti-quarks, showing that they 
both belong to the same universality class. This gives a strong argument in favor of universality 
for arbitrary products of characteristic polynomials in the regime of weak nonhermiticity . For com- 
plex models without chiral symmetry this follows from [2B] by proving universality of the kernel and 
polynomials there. 

At strong nonhermiticity we found that the partition functions of the two models can be mapped 
onto each other, provided we consider an equal number of quarks and conjugate anti-quarks. This 
correspondence indicates possible universality of some strongly nonhermitian partition functions. An- 
other argument could be that for the second model its weakly nonhermitian (universal) partition 
function can be matched to the strongly nonhermitian one by taking the limit a — > oo. However, we 
found that two partition functions with a general flavor content, for example containing only quarks, 
disagree in general. Despite the fact that they have a very similar structure, only the leading order 
term in expansion with respect to /i can be put in correspondence. This indicated that universality at 
strong nonhermiticity regime may be more subtle. When discussing such a disagreement some caution 
has to be added. As we have already remarked, a large- N limit of the partition functions with general 
flavor content as functions of the masses does not exist, properly speaking. 

Let us come back to the question of universality of correlation functions at weak nonhermiticity. 
In [5] the corresponding microscopic spectral density was computed for the first model eq. ()2.1j) at 
v = 0. There, an exact replica method was used, that expresses the density in terms of a product and 
a ratio of characteristic polynomials, 

P/| W eak(0 = ilCl^f^^^^^CD^^^CD, (3-57) 

written in terms of the rescaled variable £ = 2N(qq)z. Here by a negative number Nf = —2 of flavors 
we indicate the corresponding number of inverse powers of determinants m = n = 1 in eq. ()2.1Uj) . The 
latter objects in the present context are also frequently called the bosonic partition functions. Let us 
compare this expression to the result of [T%] 

Pn\ weak (0 



where in the second step we have used eq. (J2.18JI for n = m = 1, corresponding to eq. (|3.55j) in the 

complex plane. As we know, the two fermionic partition functions, that is corresponding to positive 

(Nt v) 

value Nf = 2, Zj jf agree for any value of v. We therefore can compare the weight function of 



= lim w(z, z) kn(z, z* ) 

lim w(z,z*) h N \Z { ^ f=1+1 ' u) (T;z,z*) , (3.58) 
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the second model in the weak nonhermiticity limit with „ 
While for the weight factor we have 



i|CI 2 ^ /= 



-1-1,0) 



that plays the same role. 



-i 



lim w(z, z*) h N 

N— >oo r— >1 



2v+l 



cxp 



(3.59) 



in terms of the rescaled variable, for the bosonic partition functions taken from [5] it follows 12 



C Zi (a;CiC ) ~ ^-exp 

or 



—(C 2 + C 2 ) 

4a 2 ^ ^ ' 



-(g)- 



(3.60) 



In the limit of small a or large \(\ 2 we can expand the if-Bessel function Kq{x) ~ e~ x /^fx and both 
expressions eqs. (j3.59|) and ([3.60[) will coincide for v = 0. However in general they disagree, leaving 
different possibilities for the universality of spectral correlation functions at weak nonhermiticity. We 
are left with two possibilities, one being that the models are in different universality classes for the 
spectral correlations, despite agreeing for all correlations of products of characteristic polynomials. 
Or, alternatively, there is no spectral universality at weak nonhermiticity at all, in the sense that such 
a universality only holds for characteristic polynomials and not for the weight function itself. 

To further illustrate these possibilities let us point out a major difference from the Hermitian large- 
N limit. At zero chemical potential the weight function corresponding exp[— NV(x)], with V(x) an 
even polynomial, is reduced to a trivial unity factor in the microscopic scaling limit Nx = const . For 
that reason the distinction between bare and full kernel becomes immaterial in that limit. Therefore 
the agreement between the two models for the product of two characteristic polynomials leading to 
the kernel, eq. (|3.55|) implies the same agreement for the spectral correlations. However, at weak 
nonhermiticity this is not the case any longer. As we see in eq. ([3.59J) the weight function remains 
different from unity in that limit. In fact it has to be a function of the imaginary part 3 : m(C) and the 
rescaled chemical potential a, such that it reduces to a 5— function in 9m(£) in the limit a — > 13 . 
Otherwise the model would not reduce to the chUE as a model of QCD pQ. 

Finally we would like to mention that there exists a matrix model [48| different from the two 
models studies in this paper. It is given in terms of two independent complex matrices and can be 
diagonalized to give a complex eigenvalue representation. The weight function of the eigenvalues of this 
model is precisely given by eq. ([3.60ft . The solution of the model |lHj by the method of orthogonal 
polynomials provides an very interesting third alternative model for QCD with chemical potential, 
having simultaneously a matrix and complex eigenvalue representation. While the microscopic density 
at weak nonhermiticity agrees with that found in |Hj the strong limit has not been studied so far. 

Another approach to the open problem of universality would be to compute the bosonic partition 



function in the second model Z 



(N f =-2,0) 
II 



Thus it remains to be seen if the universality found for 



products breaks down for inverse powers of characteristic polynomials. 



4 Phase transitions 

The aim of this section is to study the chiral phase transition at a critical value fi c . Since the second 
model eq. ()2.5|) is always in the broken phase JH]) with a constant macroscopic density on an ellipse 
for all allowed values of r £ [0, 1], we will only study the first model eq. ([2.1[) . Furthermore, we will 
restrict ourselves to the presence of quarks alone. This is mainly because only in this case a very 
compact expression eq. (j2.4[) is available (compared to eqs. (|2 . 1 1 f) or ([2.12(0 . 

12 Note that in [S| the matrices are chosen to be antihermitian. 

13 We would like to mention that in the analysis |23| only such deformations of the weight functions were studied which 
keep a Gaussian representation of the J— function. Such deformations still do enter in the macroscopic density p(0). 
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The same model has already been studied for one flavor Nf = 1 in great detail ^2]. The virtue of 
the expression eq. (|2.4|) is that most of this analysis carries over to several flavors. An extended version 
of the model eq. (|2.1j) including temperature has been studied to predict the phase diagram of QCD 
for two light flavors 9 . However, in order to be able to solve the model rather strong assumption had 
to be made on the nature of the saddle point. One of the main goals of our present analysis is to see 
if such a consideration for zero temperature can be checked and further extended, with an additional 
number of not necessarily degenerate flavors included into the model. More recently the model with 
temperature and chemical potential |Hj has been extended to study the different effect of baryon and 
isospin chemical potential for two flavors 10 . The latter corresponds to having a pair of a quark and 
its conjugate, and the authors find a doubling of the critical line as compared with 

Here there is an important difference from the previous section in the large-iV limit. We are 
interested in finding discontinuities of the partition function eq. 1)2.4(1 as a function of masses and 
chemical potential \x when iV becomes large. In contrast to the weak and strong nonhermiticity limit 
considered before, we therefore will not assume any scaling of the masses m/ with N. This will modify 
the respective saddle point eqs. ()3.5|) and ()3,31|) as the Bessel function will now make a nontrivial 
contribution to it. It introduces a mass dependence into the saddle point that makes an analysis of 
the effect of light versus heavy flavors possible. 

Before taking the saddle point limit we further compactify the expression eq. (|2.4I) . by multi- 
plying the Vandermonde determinant of the radial coordinates, Ajy / (r 2 ) = det^.j [rf^ with the 
determinant of Bessel functions. The resulting matrix elements of the single determinant read 

^2k=i r k^ * Iv{ m i r k)- Due to the symmetry of the partition function under permutations of the 
we can follow the same steps as described after eq. (|3,22[) , This reduces the determinant to Nf\ times 
a simpler determinant with elements l ^ I v (rriirj). We can now take the integration over each drj 
into the j— th column, and obtain up to the symmetry factor Nfl 

^W/}) = e-^> 2 ^-^ (4.1) 



x det 

i,j=l,...,N f 



oo 

drr ..+l + 2(i-l) (r 2 _ ^)N e -N(q q )^j^ 2N{ - q) 2 mir) 



In this form valid at finite-A?" the partition function looks exactly as a determinant over finite-iV 
partition functions of a single flavor, apart for the difference between the index of the Bessel function 
v and the different power in r to {y + 1 + 2(j — 1)). At the saddle point this difference will be of 
course subleading and allows for the analysis of discontinuities of the Nf flavors in terms of a one- 
flavor partition function elaborated in |12j . We can now evaluate the saddle points of the integrals 
individually. Taking the asymptotic limit for the Bessel function I u (x) ~ e x /y/x we obtain for each 
row 

_ r (q q f + (^)2 m . = o i= i Nf . (4.2) 

In order to lift the resulting degeneracy of the determinant we have to differentiate the Bessel functions 
as usual. Furthermore, we note that since it is only their asymptotic exponential behavior that enters 
the saddle point equation we interchange the differentiation and the saddle point procedures. The 
result is 

z\ N ?' v \n-{m f }) = 6 " ..det mf'C^f^^m^Isp , (4.3) 

A Nf (m z ) i,3=l,.. .,N f L 

where we have defined the partition function with its trivial exponential mass dependence removed, 

zj^m) = e +N ^ 2m2 Zf f=1 ' v \n;m) . (4.4) 
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The large- N partition function is thus given by a determinant of (differentiated) single flavor partition 
functions at their saddle point value. This enables us to draw some general conclusions. 

A first order phase transition occurs if the first logarithmic derivative of the partition function is 
discontinuous at some value. We look for a discontinuity with respect to so we define the quark 
number density 

^cynZf^VWI) • (4-5) 



n a 



Applying this to our result eq. (|4.3|) 

n„ = Tr In 



A 



{A)ij = m\ 1 & nii i z I (ii;mi)\ S v ) , (4.6) 



we observe that a discontinuity occurs whenever an individual differentiated matrix element becomes 
singular. Thus, the phase transition arises when the one-flavor partition function Zj 1 (fj,;rrii) with 
the smallest value [i c becomes discontinuous. Since fi c is a functions of the mass mi, we have to 
compare the m^-dependence of the saddle point equation <|4.2|) with its corresponding value fi c (rrii). 
We will find that it is always the smallest mass (which may be zero) that has the smallest value of \x c 
and thus drives the transition. 

The analysis of the saddle point solution and the corresponding partition function for Nf = 1 has 
been made already in great detail in 12 , and we follow it closely. The saddle point equation Q4.2|) is 
of third order and thus may have up to three real solution. We begin with the simplest, massless case: 

= r r {m = 0) restored 

r (qq) 2 = => r| S p = < +y/ (qq)~ 2 + fJ, 2 = r^m = 0) broken . (4.7) 



r 2_ M 2 ~ r m/ - u 'isp 



-V(qq)- 2 + v 2 i [o,«0 



It is easy to see that the solution with r| S p > belongs to the broken phase (b) with an exponentially 
suppressed partition functions, while the solution with r| S p = corresponds to the restored phase (r) 
|12j . The negative solution is rejected being outside the integration domain. If we switch on a mass 
m the signature of the saddle point will remain the same 14 , with solutions being < r r (m) < r b (m) 
as we will see below. 

The equation that determines the critical value \x c is given by the requirement of partition functions 
at two competing saddle points being equal, 

( r 2_ M 2) e < 5 g> 2 (2mr 6 -r 6 2 ) = ^2 _ r ^ e (q q ) 2 (2mrr-r?) _ ^ 

In the massless case m = this leads, after inserting the solutions (|4.7|) . to the following equation for 
the critical line: 

i + (qq) Vc + H(qq) 2 ^} = o , => (qq)» c ~ o, 527 . . . , (4.9) 

where we have given the approximate numerical value for its real solution (see also fig ^ below). 

Let us now determine how this value shifts if we include a small mass. We first present a pertur- 
bative analysis to the leading order in the mass 15 and then come back to the full solution below. First 
we determine the shift of the saddle point solution due to the mas term, 

r r (m)-r r (0) = 5 r {m) - m ^)V 



1 + <gg> V 



Tfl 

r b (m) - r 6 (0) = S b {m) = - ^- g p 2 ) ' ( 41 °) 



24 



Hc{m) 



"TO 

m 



Figure 1: The critical potential as a function of mass, both given in units of (qq). The critical value 
at m = is /2 C (0) = 0.527. . . . 

where we retained only the leading order linear behavior in m. 

Next we compute the resulting shift in the critical line, determining fi c (m): 



1 + (m) Vc + H(qq) 2 £] - 2m(gg)V<?9)- 2 + ^ + <?(™ 2 ) = . (4.11) 
To the leading linear order in m both S r ^ dropped out. We thus obtain for the shift in ^ defined 

as 

£{m) - £(0) = tH , (4.12) 

the positive quantity 

= 2(qg)m(qq)^ c (0) 2 
VI + (??)Vc(0) 

after inserting it into eq. (|4.11|) . It is again given to leading order only. From this result we can 
deduce the following. If we have massless and massive flavors present in the partition function eq. 
(|4.H|) the critical value for the massless single flavor partition functions will be reached before that of 
the massive flavors, as /U c (0) < \x c (m ^ 0). Second, if only massive flavors are present it is the critical 
Hc{m) of the lightest flavor which is reached first, as for m\ < m2 holds ^(mi) < /^(m^). In both 
cases it is the lightest or zero mass that triggers the phase transition. Since it is of first order for a 
single flavor ^] it is thus first order for any combination of massive and massless flavors, as follows 
from eq. (|4.6|) . Moreover, for each different mass there is separate transition, leading in principle to 
a sequence of transitions 16 . A similar feature of two distinct first order lines was found in ^01 for two 
flavors. 

We have also convinced ourselves that the monotonic behavior of the critical chemical potential 
/U c (m) persists beyond the linear approximation, eqs. Ij4.ll-} (I and (j4.12j) . This can be done by inserting 
the two positive real solutions for the third order saddle point equation into the condition eq. (j4.8j) 
and solving it numerically. If we measure all quantities in units of (qq), f = {qq)r, fi = {qq)fJ- and 
= (qq) m ( or set {qq) = 1) we have 



m 



-V^l + Zi 2 ) +m 2 cos[<^/3 + 4tt/3] + rj 
o o 



14 This can be seen from the discriminant for all real values of the masses. 

15 A similar analysis was made in for Nf — 2 with a different chemical potential for each flavor. 

16 We wish to emphasize that in general the flavor dependence of matrix models is too weak. In QCD the sign of the 
/3— function changes for sufficiently many flavors, and the finite temperature transition changes from second to first order 
when changing from 2 to 3 flavors, which is not reproduced from the matrix model U3- For that reason our conclusion 
should be taken seriously for small Nf only. 
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Figure 2: The logarithm of the partition function for N=10 (top), N=20 (middle), and N=30 (lower 
line). The vertical lines indicate the positions of the of the corresponding critical values ^(mi^). 



m 



n 



-V3(l + /i 2 )+m 2 cos[v9/3] + - 



COS [if] 



m(m 2 — 9/j, 2 + 



(4.14) 



[3(1 + ft 2 ) +m 2 ]2 

as the two positive solutions of the saddle point equation. At m = (<p = tt/2) the first solution f r 
vanishes as it should. We have plotted above the numerical solution of eq. (|4.8|) . after inserting into 
it the full solution eq. (|4.14l) . It can be seen that /U c (m) is a monotonous function beyond the regime 
of small masses m. 

In order to check our findings we have examined numerically the partition function for two flavors 
Nf = 2 with two nondegenerate masses, in the representation of eq. (|4.1|) , We have chosen the 
parameters mi = 1 and m<i = 2 in units of (qq) and v = 0. The resulting logarithm of the partition 
function is plotted above for three different values of N as a function of [i (figure • Even for the 
smallest value of ./V the points where the derivative of the partition function and thus the free energy 
becomes discontinuous are very well visible. The positions coincide with the critical values of \x for 
each mass scale, p, c (fh\) = 1.2119 . . . and fi c (jn2) = 2.119 .... (from fig. Q). 

After having made a general statement we give an explicit example for the partition function eq. 
(|4.3j) below and above the transition. Here we can again make use of analysis |12j where we briefly 
repeat some of the results. The single flavor partition function (modulo its mass prefactor eq. (|4.4l0 
can be conveniently written as 



zi(\x\m) 



dr r(r 2 



fi 2 ) N e- Nr2 I (2Nmr) 



+ 



dr r(r 2 



fi 2 ) N e- Nr2 I (2Nmr) 



z r (fi;m) + z b (fi;m) . 



(4.15) 



We used |12j that the restored symmetry solution satisfies r r G [0, // 2 ] and the broken symmetry 
solution is r b £ [fj, 2 , oo) (see also eq. H4.14(l ). Here and in the following we set (qq) = 1 and v = as 
it is done in ^2]- The two different solutions can be written as a double sum at finite- N, 



z r (n; m) 
z b (n;m) 



2(^+1) „ M 2 



{N + s)\ 



(Nm 2 ) l (Nn 



l\s\(N + I + s + 1)! 



EE 

1=0 s=0 



,2\l+s 



1=0 s=0 



l\s\(l-s)\ 



(4.16) 
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These formulas easily generate the massless partition function and its derivatives, z^ b (fj,;m = 0), as 
we need them in eq. (|4.3|) . We note that only even derivatives contribute and give the examples 
necessary for Nf = 2 massless flavors. They read in the restored phase 

z r (fi;m = 0) ~ /i 2(JV+1) 

4'fo™)L=o ~ M 2(iV+1) 2iV M 2 , (4.17) 

and in the broken phase 

z b (n;m = 0) ~ e"^ 2 

4V;™)U=o ~ e-^ 2 2iY 2 (l + / u 2 + l/iV) . (4.18) 

Taking the limit of two degenerate, massless flavors in eq. (|4.3j) we obtain for large N 

< N =2 ,v=o), » / 2/(A«;0) \ / n 4N restored 

^ ^ m = 0) = det [ zJ'( M ;0) J~\ e- 2 ^ 2 (l + / , 2 ) broken " (419) 

This compares to the single flavor result 

(jv / =i )l /=o) / . n x / ^ 2N restored 

^ (M;m = 0) ~ | ^ brQken . (4.20) 

We see that the behavior is very similar, with the quark number density in the two flavor case being 
twice as high. 

As we have mentioned already there is more than one transition in the presence of several mass 
scales. In particular this opens the possibility to have for example two transitions for two massless 
flavors and one massive flavor, as it is often used in a simple model for QCD. For /U c (0) < /x < /j, c (m), 
the massless building blocks are already in the restored phase while the massive ones are still in the 
broken phase. This is an intermediate regime before full symmetry restoration. It would be very 
interesting to use the exact solution of our model as a further testing ground for lattice algorithms 
as for example in F° r the discussion of phenomenological consequences we refer to [10] where a 
similar phenomenon is observed. 



5 Summary 

We have computed and compared two different matrix model partition functions for QCD with chemi- 
cal potential. While the first model has only a matrix representation and no eigenvalue representation 
the situation is the reverse for the second model. We give very compact, new expressions for finite- N 
for the first model in terms of integrals over the radial coordinates of the underlying sigma-model-like 
representation. They hold for an arbitrary number of quarks or both quarks and conjugate anti-quarks. 
By taking the large-iV limit these expressions are then compared to the corresponding results from 
the second model expressed in terms of orthogonal polynomials in the complex plane and their kernel. 

In the limit of weak nonhermiticity we find a complete agreement, with a different identification 
of the nonhermiticity parameters fi and r and mass parameters in the cases of only quarks present, or 
both quarks and conjugate anti-quarks. 

The matching between the two models in the limit of strong nonhermiticity holds only for an equal 
number of quarks and conjugate anti-quarks. For such a matching the masses have to be rescaled 
with respect to the \i- or r-dependent macroscopic density proportional to the level spacing. For a 
generic number of quarks and conjugate anti-quarks, including the case of only quarks, a mapping 
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can be achieved when perturbatively expanding in small chemical potential //, despite the very similar 
structure of the two partition functions. However, in this case a proper large-iV limit of the two 
partition functions does not exist. 

From these finding we concluded that at weak nonhermiticity and also in a special case at strong 
nonhermiticity the expressions for the fermionic partition functions are universal. The issue of uni- 
versality of eigenvalue correlation functions is left open, as it involves also the computation of bosonic 
partition functions. 

In the last part we have investigated the phase structure of the first matrix model, which is known 
to possess a first order phase transition for a single or several degenerate flavours (the second model is 
always in the symmetry broken phase by construction). We have analyzed the model for an arbitrary 
number of nondegenerate quark flavours, exploiting the compact expressions derived for finite-iV. 
In the large- iV limit they give the multi-flavour partition function as a determinant of single flavor 
partition functions of different masses, and their derivatives. Consequently for each nondegenerate 
flavor a separate first order transition occurs, where the critical fi increases monotonically in mass. 
We have checked our findings numerically for two flavors with nonvanishing, nondegenerate masses. 

It would be very interesting to compute also the corresponding bosonic partition functions and 
compare them, in particular in the light of their relation to complex eigenvalue correlations and the 
related question of universality. 
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A The sigma-model representation of the partition function 

Here we derive the sigma-model representation for both the partition functions eqs. (|2.5j) and (|2.1[) , 
The latter can be obtained by omitting all conjugate flavors. To this aim we first unify notation, by 
writing eq. (|2.10f) as 



as well as choosing signs for the chemical potentials /if = +fi for / = 1, ... ,m, and fj,f = —jj, for 
/ = m + 1, . . . , Nf = m + n. After introducing two sets of complex Grassmann vectors \A and xb of 
size N and N + v respectively, we can express the determinant in eq. (|2.1|) as 




(A.l) 



denoting the complex masses with 





"l 



(N f ,u) 



/ 



d 2 XAd 2 XB e 



£ flJufiXA jXAf+Xa fXB f) + m/(x a/ 1 




') 




(A.3) 
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We define the expectation value over the ensemble amounting to the integration over matrices as 



It can be further rewritten as 



iE/ii(xi f®XBf+X B f&XAf) 



exp 



exp 



Tr(*E/iiXs/®x^/ + ^Y^UxAf^Bf) 



N r 



iV(OT} 5 



Y, ^(XA f ®X ] Bf ){XBg®x\ g ) 



/,<?=! 
Nf 



Yl (XAgXAf)(XBfXBg) 



/,<?=! 



dQclQ^e 



\jf f ! g= x{Sx\fXA B )Q^+{x^ Bf XB g ){Q^ f ) -N{qq)*TTQ\Q 



(A.5) 



where we have integrated out the matrices <£, performed the trace and made a Hubbard-Stratonovich 
transformation. The auxiliary matrices Q are now quadratic of size Nf x Nf. The integrations over the 
Grassmann vectors \A and xb in eqs. (|A.3|) and (|A.5|) are now all Gaussian, leading to the following 
determinant 



dQdQ ] det 



(M + Q) T ®1 N fiS 3 ®l 



N 



(M + Qt)T g, ljv+t 



-jV(g 9 } 2 Tr QtQ 



dQdQ 1 " det [M + Q f ] 1 



(A.6) 



x det 



(Af + Qt)(M + Q) - /x 2 (M + gt)S 3 (M + Q^y 1 ^ 



N 



-N{qq} 2 TrQlQ 



Here we have introduced the mass matrix M = diag(tJi, . . . ,u>N f ) and the generalized Pauli matrix of 
size Nf x Nf 

' lm 

-l r 



E, = 



(A.7) 



Furthermore we made use of the following property of determinants for invertible matrices D 



det 



A B 
C D 



det[£>] det [A - -B D~ C\ . 



(Ai 



For real valued masses, M = Aft, we can further simplify eq. I)A.6|) by shifting Q + M — > Q, 



-JV<OT> 2 Tr(Q+Q - M(Q+Qt)) 



V 



xe 



(A.9) 

In the absence of conjugate anti-quarks, S3 = 1jv>j the second determinant further simplifies, leaving 
H 2 lN f only. In order to arrive at a unitary group integral we introduce the "polar decomposition" 
Q = U R, where W = U~ l is an Nf x Nf unitary matrix and W = R is hermitian and positive. 
The Jacobian for the transformation is derived in the following appendix. An alternative is the Schur 
decomposition of the matrix Q used in which quickly becomes cumbersome for more than Nf = 3 
flavors. The advantage of our present method is that we can conveniently exploit the saddle point 
approximation for the integrals over positive eigenvalues rj of the matrix R. 
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B Jacobian for the polar decomposition 

A classical result of linear algebra states that any Nf x Nf complex matrix Q can be written as 
Q = UR where = U' 1 is an Nf x Nf unitary matrix and W = R is Nf x Nf hermitian positive 
matrix. That is the matrix-generalization of the "polar decomposition" of a complex number q = re id 
Clearly enough the number of real degrees of freedom of the matrix Q (i.e. 2Nj) matches with the 
total number of real degrees of freedom of U (i.e. Nj) and R (i.e. Nf). Under this matrix change of 
variables, the integration measure produces a Jacobian: 



dQ dQ^ = J(R) d^i{U) dR . 



(B.l) 



where dfi(U) is the Haar measure on the unitary group, dR is the measure on the space of Hermitian 
positive matrices and J(R) is the Jacobian that we are going to calculate in this appendix. For 
calculating the Jacobian 17 , we differentiate: 



dQ = dU R + UdR. 



(B.2) 



As R is an Hermitian matrix, it can be diagonalized by an unitary matrix V: R = VfV~ l with f = 
diag(ri, . . . ,r n ) having positive entries. If we left-multiply eq. 1)B.2|) by and right-multiply 

by V, we get 

V~ l U~ l dQV = V- l U~ l dUVr + V~ l dRV . (B.3) 

This equation can be written as dQ = dU r + dR with dQ = V^XJ^ 1 dQ V, dU = V^ 1 !/^ 1 dUV and 
dR = V- l dRV. The matrix dU is anti-hermitian, since (cfl7) f = V ] (dU)^UV = -V^U^dUV = 
—dU (we used dU^U = —U~ l dU). Moreover, the matrix dR is hermitian, since dw = (V" 1 dRV)^ = 
R. Therefore, the relation between the differentials of the independent complex variables is: 



(dQ)u 
(dQ)ij 
(dQ)ji 



dUu Ti + dRu 
■ dUij Tj + dRij 



Let us introduce a real parameterization as follows: 



l,...,N f 



i < j 



(B.4) 



(dQ)ij 
{dil)i 3 



dxij + idyij , (dR)i 



{dR)ij 
(dR)u 



dpij + idqij 
dpu 



(dU)ij 
(dU)u 



dvij + idwij 
idwn 



1 < J , 

i = l,...,N f 



i < J , 

i = l,...,N f 



By matching the real and imaginary parts we have: 



dxi 
dx; 



dpu and dya = ndwa, for i = 1, . . . , Nf. The corresponding Jacobian is J\ 



n2W 



is 



dvijrj + dp^ and dxji 



-dvifTi + dpij for any pair i < j. The corresponding Jacobian 



det 



d(vij ,Pij) 



det 



3) dyij = dwijVj + dqij and dyji 



Js = IIS det 



-1 



n 



dwijTi 
i<j{ r i+ r j) 



(rj + rj), and therefore J 2 = Y[f<j( r i + r j)- 



dqij for any pair i < j. The corresponding Jacobian is 



7 For general methods of evaluating Jacobians, see |49| |5(J|. 
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Therefore the full Jacobian is J = Y[f=i r i Ylf<j( r i + r j) 2 - Going back to the original variables 
dU — > (iii — > <ii? and dQ — > <iQ does not produce any additional factor in the final Jacobian (up 
to an overall sign) as they differ only by unitary transformations. We thus obtain the final result 

N f Nf 

dQdQ^ = YlnYlin + rjfdfx^dR 

i=l i<j 

N f N f Nf N f 

« II n II ^ + r tf 11^ - r rf MUW(V) n dn 

i=l i<j i<j i 

Nf Nf 

= \[r i dr i \[ = {r 2 i -r1fd l x{U)d l x{V) , (B.5) 

i=l i<j 

where we used the standard fact that R = VfV~ x implies dR oc Y\i<j( r i ~ r j) 2 Y[f f dr{dfj,(V). We 
have thus shown eq. (|B.5|) to be the expression for the Jacobian for the polar decomposition Q = U R. 
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